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Symbols 



E 5 = Vector magnitude of incident field 

B g = Geodesic field vector 
)Co= to/c = Free space wave number 
^ - Relative permeability constant 
c '= Complex relative dielectric constant 

c '= Complex relative dielectric tensor 



n = Complex vector index of refraction 

y = j kb n = Complex vector propagation constant 
I 0= Angle between yand Bg 
? - Cos i j 

X/ V/ Z = Cartesian coordinates 
f/ 0,<p= Spherical coordinates 

Small displacement along X-axis 

F = Force vector 

e = Positive magnitude of electron charge 
m = Mass of electron 
- e/m = -1.77 X 10^^ coul/Kg 

i;= Average number of electron- ion collisions per second 
r = Radius vector 



V = Vector velocity of an electron 
N = Density of electrons per cubic meter 

_7 

jig = 4 77 X 10 h/m = Permeability of free space 
1 ~9 

fo *=« X 10 fd/m = Permitivity of free space 
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= Resonant angular frequency, lossless plasma 



eBg = Cyclotron, or gyro, frequency 
m 

E, H = Electromagnetic field components 
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Symbols (Continued) 



D = Electric flux density vector 

t = Current density in amps, per sq. m. 

o = Surface charge density in coul. per sq. m. 

• s 

p = Charge density in coul. per cubic m. 

X = t0p2/ ^2 

Y = ^/ + Y ^ 

X y z 

Y = Y cos Q, Y = Y sin 6 cos <o, Y = Y = Y sin 0 sin <a= Y cos tb 

X y z L ^ 

= 0' + ";2' =*0TT2 =Ysln# 

X Y L 

r = i-j i^/co 



= n^ - 1 = n . n -1 



a . = Unit vectors 
1 



A,=a^xa =a, sin r 
1 6 z 1 ^ 



X = 1 - r X 




Z= 1 - X 



E 

xy 



= '^ 777 ^ 



* 





E 

X 



= . Critical frequency 

to = Ordinary wave critical frequency 

to , = Left circularly polarized wave critical frequency 
cL 



to „ - Right circularly polarized wave critical frequency 
CK 

o = Lorentz conductivity tensor 



p = Lorentz resistivity tensor 
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Symbols (Continued) 



I = Identity tensor 

s = Length of path through the ionosphere 
0 

c « 3 X 10 m/sec = Speed of light in a vacuum 

V = Group speed 

g 

V = Phase speed 
P 

Oi = Attenuation constant 

j3 = Phase constant 

do = Angle of incidence 

^ = Angle of refraction 

R = Measure of polarization 

n.. = n, n. = Components of refractive tensor 
1] i ] 
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1. PIASMA FREQUENCY 
1.1 Lossless Case 

Consider a lossless plasma having a density N 



vt^ 



of free electrons. If the electrons are compressed 
_^X laterally through a small displacement, 4 / along the 
X-axis, a surface charge density 



x=o X= C 
Figure 1 



-P 3 = -Ne 4 



( 1 ) 



may be considered to exist at y = 4 • Then, from the definition of the 
electric displacement vector, from y= o to y= 4 / 

P<=* 



Hence, 



Dx=^oE^=NeC ' 



or the force F is given by. 



F = -eE = - Ne 

X X “7" 

^0 



c* 



Since force is given by the product of mass with accelleration. 



m d 4 
dt' 



+ Ne' 



^ = 0 . 



( 2 ) 

(3) 

(4) 

(5) 



2 2 

Substituting o) = Ne into (5), the harmonic equation for free oscill- 

m f 0 

ations of the electrons is obtained, iTamely, 






dt 



T 



+ tOp 4=0. 



( 6 ) 



The solution of (6) is sinusoidal with an angular frequency co . Thus to 

P P 



is called the resonant angular frequency for the lossless plasma. 



1. 2 Lossy Case 

Suppose there are an average number, u , of collisions per second 
between electrons and positive ions. This introduces an average loss of 
momentum, or an equivalent fluid resistance term, into equation (5), which 
equation becomes , 

2 2 
m d ^ + pm d 4 + Ne t = 0 . (7) 

dt 2 dt c 0 

pt 

Dividing by m, and using the operational method of substituting 4=4o® 



1 



4 





1 




J 



m 



1 



into equation (7), the determinantal equation is obtained, 

(p2 + i/p + to 4 = 0. (8) 

Equation (8) yields the characteristic roots , , 

p = -1//2 ±j tOp V 1 - ( (9) 

From equation (9), it may be concluded that ionic collisions in the 
plasma introduce a time decay constant of 2/p, and reduce the plasma 
frequency by a factor V ^ ^^^2 . For a low loss plasma, the 

latter factor usually is ignored. 



2. HOMOGENEOUS ISOTROPIC IONOSPHERE 
2.1 Low Loss Case 

Consider a small collision frequency for a homo- 




geneous isotropic ionosphere of electron density 

N(2i_) . Let a uniform plane wave front, Eope^^^^ 
m3 

be incident upon the ionosphere at an angle of 



incidence 0^. If r is a radius vector describing the 
path of motion of an electron having a velocity V, the force F upon an 
electron is given by. 



F = -eE = 91 



d.2 



+ p m d r 
dt 



or 



d V + p V = -_e_ E 
dt m 



Substituting V e into (11) and solving for the steady state velocity. 



V = - e E 



= - e E 



( 10 ) 

( 11 ) 



f ^ \ • / • 1 • ( 12 ) 

m (joH- V) jm( to -] V) 

The equivalent current density due to the motion of electrons then becomes, 

7 = p V = -NeV = Ne^ E 



/ . X • (13) 

m (p+ jto) 

Equations (12) and (13) show that the electron path, velocity, and current 
density are linear in the direction of E. 

In order to see the effects upon the electromagnetic wave. Maxwell's 
equations in the steady state phasor form will be used. 



2 



I 




The Maxwell's equations are, div B = 0, div D = 0, curl E = - b / 

Bt 



curl H = t / B = /Li' jU 0 H , D = f ' fo E . 

Bt 

The corresponding complex form is, 

V. H= 0 , V.E = 0 , vxE = -ja)jUoH, 
V X H = L + jo 5 €o E. 

Substituting from equation (13), 

2 r 



(14) 



(15) 



vx H = Ne E + joJ CqE = jCL'€o A" 
m (p + jto) 

[1 + ^ 



(16) 



or. 



joJCo m(y+ jio) 



VxH = j 03 € o[l + 



0 ), 



2 ii>- jitj) 



) 






jco €, 



( 1 - 



COr 






+ 03^ 



U 60 ^,^ 1 - 

-3 7^ ^ E , (17) 

+ 03 ^) J 



2 2 

in which to = Ne was also substituted in the right member of equation 

^ m c o 

(17). 

obta ined , 



From equation (17), an equivalent complex dielectric constant € is 

c 



U3r 



cl =€ '-3 €" = ( 1- 

c 2^ 2 

03 +P 



_i P 03 



J -3 



/,2 . 2 
03 + y 



(18) 



Returning to euqations (15), the wave equation will be determined 
by elimination. 

- - 2 “ 2 - 

vx ( vx E) =V( V.E)-V E = -v E = -jo 3 /Uo*' 



V X H . 



(19) 



Using y = j j3 = jk^ V c ' , and = 03 V jo o € 0 = 0 )/c , 

c 

substituting for vx H, equation (19) becomes, 

v^ E + k^ € ' E = 0. 

° c 



upon 



( 20 ) 



For V «o, n ='\TP will be evaluated. 

n = (e'“3€^) = '/ 



: V e’ ( 1-j -S — ) = V c 
2 c' 



,, 1/2 

( 1-3 _£_) 

_ €' 

'’trT~ 



(21) 
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and hence, 

y = j k„n « + j ko7"r~ . 

2^/7T 

From equation (18) with 



ypa V ko 



+ Jk„Vl-(^)2 

^ 2 W 



to 



The phase velocity is. 



V s to = to/ k. 
P 



^ V 1- ( ^p)2 V 1- ( t0p)2 



to 



to 



and the group velocity is , 



to 



( 22 ) 



(23) 



(24) 



(25) 



2.2 Lossy Cas^ 

For the losiy case, the exact value of n = Ji* should be used. For 
this purpose, consider \/'a ± j ^ ' , B = tan b/a . 



a- 



Va ±j b = 1/a2 + b^ e 

‘ Vi: 



cos 



a“^ + b^ 



sin B/2 = 



7? Va2 + b2 . 



Va ± i b = L- V 






+ a 






'rr 



+ b^ - a 



Substituting from equation (18) into (26), 



1 



r 



to 



+ ( 1 - 



2 + 2 

to 



-d- V 



2 ^ 2 

to + p 



1/2 1 

^/r. 

1/2 



(1 



- 



2 + 1/ 



2 tOp"^ 

) + 



60 



(1- P 



2. 2, 2.2 

to ( to + P ) 

2 to 4 ^2 
) + P 



(26) 

1/2 

1/2 



0? + p^ 



2/ 2 ^ 2^2 
to ( to + p^) 



(27) 
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Thus , y = ko [ n]and /3 = k^Re [ n]. (28) 

The Re [ n]is also spoken of as the refractive index for the ordinary wave. 



2.21 Critical Frequency 

From equation (28), the phase velocity within the ionosphere is. 



P ./ / — 

^ € 0 Re [n] (29) 

Designating the angle of refraction by 0^. , from Snell's law, 

sin 6 = sin ^ ° 

^ Re [ n] (30) 

If 0r < 1 t /2, the wave enters the ionosphere, is refracted, and passes 
through the ionosphere, unless the right member of (30) is greater than 




If the right member of (30) is greater than unity, no real angle 0^ exists and 

the wave is reflected. Designating the critical angular frequency by o) , 

c 

this frequency is determined by 0^ = n/2. Thus, for any angle of incidence, 
the critical frequency is determined from. 



s in 00 = Re f n ] . 


(31) 


For the low loss case, from equations (23) and (31), 




sin e, = Vj. ( Wp ,2 _ 


(32) 


CO 




c 




05 = CO sec 0 o 

‘"n n ° 


(33) 



Equation (33) is known as the secant law, and yields w = co for 

c p 

normal incidence. Substituting the numerical values of e, m, and into 
(33), 

f = 9 ^/iT sec 0 0 , (34) 

c 

in which f is in cycles per second. If the electron density is given per 
c 

cubic centimeter, the f will be in kilocycles per second. 

c 
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At a given angle of incidence, frequencies f > f will pass through 

c 

the ionosphere, and frequencies f < f will be reflected. 

c 

The preceding analysis has ignored the earth's geodesic field. All 
waves acting accordingly are termed ordinary waves. 

3 . HOMOGENEOUS ANISOTROPIC IONOSPHERE 
3 . 0 Direction Cosines 

In order to take into consideration the effects of the geodesic field 
upon an electromagnetic wave propagating in the ionosphere, various 
vector directions must be taken into consideration. Hence, a table of 
direction cosines is necessary. Such a table may be formulated either 
by projections or by spherical trigonometry, or by a combination of 
both methods. Referring to the figure: 

cos 0 = cos 0 cos 7 j /2 + sin 0 x 
sin t/2 cos ( it/ 2 -(p) 

cos il)= sin 0 sin ^ 

a^ . a^ = cos i/f sin 0sin <p 

“ a , . a = sin ib 
ip z ^ 

- a , . a = -cos lb s in sf 

!/) y y o 

- a , .a = -cos )l)cos $' 

0 X 

cos 0= cos lb cos n/2 + sin lb X 
sin ff/2 cos s' 

cos Q' = cos 0 

s in 0 

sin q'=Ji-92^JT 
' sin lb 

V l-sin^Q sin 2 (q-cos2 (q - 
' sin^ lb 
sin 0COS (p 

Figure 4 sin ib 

sin ib= Vl- sin2 0s in^ (p 
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From the above, the following table may be compiled: 







a 

y 


a 

z 


■^0 


ir 


cos 0 


sin 0COS <p 


sin 6 s in (p 


0 


a0 


-sin 0 


cos 0COS (p 


cos 0 sin <p 


cos dsln CD 
sin 0 


a 

<P 


0 


-sin cp 


cos <p 


COS cp 

sin i/) 


0 


-sin6cos 0sin <p 


-sin-^ 0sin(/3Cos(/5 


sin 0 


1 


sin ^ 


sin 0 



(35) 



3 . 1 Electron Current Density 

Assume the Z-axis of a coordinate system to lie along the geodesic 

field in the ionosphere, and assume an electromagnetic field, E = E^ a 0 , 

to be incident at angle 0^ upon the ionosphere. Equation (13) gives the 

current density, due to the presence of the electromagnetic field along 

with the free electrons. However, it ignores the presence of the geodesic 

field, the geodesic field being a^ . 

The zero order component, of equation (13), reacts with to 

produce a new component / of the current density, normal to both 

Lo and B . This component lies within the XY - plane. The i , component, 
Z 1 

in turn, produces another component, u , normal to both i, and B , and 

I 1 z 

hence also lies within the XY - plane, but it is rotated t/ 2 from tj. The 
new component, , produces still another component, , normal to 
within the XY - plane. Since h^s been rotated it from f, , it is in 

the same direction as , but of the opposite sense. Also, the component 
, produced by / and is also of the opposite sense. Continuing, an 
infinite alternating series is obtained for the two directions lying within 
the XY - plane. These series will converge under certain conditions. 

They will be obtained analytically in the following analysis. 

In the presence of the geodesic field, the force equation (11), for each 
electron, must be modified in accordance with ampere's force law to 
become , 
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( 36 ) 



m y +umV+eVxB = -eE . 
dt 

Equation (36) is to be solved by iteration. That is, the solution of 
equation (12) will be substituted into e V x B and a correction term , 
obtained. The term will then be substituted into e V x B to obtain another 

correction term , etc. Before substitution, the successive cross products 

will be formulated. From table (35), 




But from table (35), 



a a sin 0 + a cos 0cos (o+ a cbs 9 
0 X y ^ z 



sin <p, and 

A, = a X a = a, sin T / C = cos 
1 6 z 1 > ^ 



(37) 



-1 






X y z 

- sin 0 cos 0COS <p cos 0 sin 
0 0 1 



a cos ©cos <n + a sin 0 

X ^ y 



(38) 



Ao= ( a ) X a = a. sin f = 

9 ' O *7 !7 / ^ 



2 ' 0 z' z 

a a a 

X y 2 

cos 0 cos cpsin 9 0 

0 0 1 



= a s in 9 -a 

X y 



cos 0 cos 



(39) 



Hence , 
and 



E = -Eo sin e , E = E^ cos ©cos (p, E = 
X y ° ^ z 

Eq cos 0 sin (p . 

A, = a E ~ a E 
1 X y y X 



A- = - (a E + a E ). 

2 X X y y 



(40) 

(41) 

(42) 



From equation (38) and (39), it may be seen that both A^ and A^ lie 
within the XY - plane, and that their slopes are negative reciprocals. That 

is to say, A^ is rotated ir/2 from Aj^. Furthermore, since each unit vector 

a. , j = 1, 2 , , is within the XY - plane , each a , , , = a . x a also is 

4 . 1 + 1 1 2 
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within the XY - plane, and a. , is in space quadrature with a.. 

1+1 1 

Incidentally, the unit vectors a , a and a . are also coplanar, with 

Z Z u 

a„ at an angle i/2 + ^with a . 

For the iteration, from equations (12) and (10), 

V = .. a ft p = e VgX B = jm (o) (43) 

° j m ( o)- ji^) 1 z I 

Thus, upon defining the gyro frequency, 

(44) 

m 

Vj= X = -e Eo 

[jm ( ui-ju) f m { ui - jv)^ 

co^E. sin C . ( 45 ) 

m ( to- j v)^ 

Repeating the procedure for 





fa 


= -e V, X B = j m 
1 z 


(to- iv ) V 2 


/ 


(46) 




^2 


= e E, B^ tOj^ Aj X 


a 

z = 










j m^ ( c»>^jt')3 










-j 


2 

e tOj^ Eg sin ^ 






(47) 






' ( to - jl^)^ 


2 






For 


^3 


= -e V„ X B = j m 
^ z 


( to- jv) V 

0 


i 


(48) 




^3 


= e t0j^3 E^sin ^ 


/ a = ♦ 




(49) 






m ( to- iv)^ 


3 3 




Continuing , 




= J e tOj^ E„ sin C 


*/i = 


- a„ , 


(50) 






m ( to - ju)^ 


4 , 4 


2 






^5 


= -e tOj^ Eo sin ^ 


a , a = 


®i ' 


(51) 






m ( to- jl^) ^ 


5 5 






^6 


6 

= - j e tOj^ Eo sin ^ 

7 


- ^6 ' ^6 = 


^ 2 ' 


(52) 



m ( to- jv) 

The above procedure may be continued indefinitely. 
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Upon summing the corrective components (43) through (52) ad infinitum. 



factoring - j co e^, and substituting o) , the velocity vector for an electron 

P 



u Ne 

becomes , 

V = + j 






+ j a, 



Ne 



6 w ( 03- ]» 



03 



2 5 

3 % 



2 

CO ^ av 



CO ^ CO ^ 
b + P o 



03 ( 03- jP)2 (ji( OD- iU)^ 






03 ( 03- j 

03 2 03 6 

p b 

03 ( 03 - jP)^ 






sin ^ ~ a. 



sin n 



03 2 03 2 

P b 



2 4 

P b 



03 ( 03- jp) 3 03 ( 03 - jp) ^ 



(53) 



or, - 

V = + j 03 € .E. J a 
Ne 



co^ 



1 



+ j a 



1 



03^4 
+ b 



( 03 - JP) 



03 4 
b 



+ 






03 ( 0^ 

sin C “ ^2 ^P 



03 ^ 03_ 



03 ( 03- 



03, 



1 + 



( 03 - j p) 



Sin ^ 



7 . 



03 ( 03 - jp)2 



1+ V 



( 03 - j P ) 



(54) 



( 03- jp)"^ 

The infinite series within the brackets of equation (54) converges for 

_ < 1 , (55) 



03b2 



2 ^ 2 
03 + P 



that is, for 



03 



> 03, '/ 1- 



( l/ 03j^)' 



Hence, for the frequency range of convergence, the brackets become. 



1 + + — 
( 03- j p)2 



= (a rjp) ^ 



^ - ,...2 ..2 . 



(56) 



(57) 



1 - 03 2 /( 03- jp) (^" jp) ^ W 



and hence. 



V = + j 03€ .E„ 
Ne 



^ 03u 03 sin f 

jg + b ^p ^ 

( 03- jp) 03[( or jp)*^ - 0^2] 






03. 



03 - jp 



(58) 
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3 . 2 Path of An Electron 



For an examination of the path followed by an electron, it is preferable 
to eliminate in equation (58), and integrate in time to obtain a radius 
vector, r , whose terminus describes the path. Accordingly, 



r = E„ a 



m 



9 



1 



+ ^b J cos^ 9 sin ^ <p 



Ui- iv) 2 2 , 

0) [ ( CO- ] t^) - CO]-, ] 






CO, 



CO- j u 

in which, from equation (37), sin ^ has been substituted in the form 

= 



(59) 



sin^ 



2 2 
cos 9 sin <p 



(60) 



Assuming i/ « 0, equation (59) shows that as the electron attempts to 
vibrate in a path parallel with the exciting field, it has an additional 
elliptical component within a plane normal to the geodesic field, and 
resonates at co = ct|^ . For frequencies co< cOj^ , the ionosphere acts as 
a conducting medium, for the current density produced by the magnetic 
field, rather than as a dielectric medium. 

The electron path, as it spirals about the direction of the exciting 
electric field, is somewhat like a trochoidal epicycloid. The electron 
current density, t = - NeV, may be considered as a source of radiation, 
a vector potential formulated, and the resulting fields computed there- 
from in accordance with Huygen's principle. Thus, the electric fields 
are no longer confined to the path of the driving field. It will be shown 
later that D, H, and yconstitute an orthogonal system, with E not 
necessarily orthogonal to y . 



3 . 3 The Curl H Equation 

The current density will be formulated from equation (58), by mul- 
tiplication with - Ne, and substituted into ctirl H in order to formulate 
the dielectric tensor. The dielectric tensor subsequently will be used 
to determine the indices of refraction. The resulting curl equation is. 
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^ X H = j c 0 CoE->( a. 



*2 



1 60 2 

1 - 



60 (cO - jv) 



- j a. 



“b “d 

2 2 

CO [ (co- ji^ ) co-^ ] 



60(C0- iu)l ( 60- iv)^ - oi^'] 



1 



( 61 ) 



The wave corresponding to the first term of the right member of equation 
(61) is sometimes called the ordinary wave, and the wave corresponding to 
the second and third terms is companionably called the extraordinary wave. 



3.31 Cartesian Coordinates 

Using E-=a E +a E +a E along with equations (41) and (42), 
Q X X Y y 

equation (61) may be re-written. 



i V X H = a 

j OJ€o X 



F 

- E b p 
X 




1 60 

1 - “^P 



60 ( 60- jv) 



- IE ^b ^p 

y 



60 jj 60 - jv)^ 



iti( iC - iv) [(60-jvf - 60 j^] 



+ a 



2. - y 



60 

1- P 



60 ( 60 - ]U ) 



. -o 60, 60 2 

J E b P 

X 



- E. 



9 2 

60, 60 
b p 



60[(60- ] 



60 (60 - jp) [(60- jp) 



+ a E 
z z 



- 1- ] 

60 (60- jv)j 



Now let 



X = 1 - 



to. 



2 2 . 
60, 60 
b P 



2 2 

60 ( 60- jl^) 60 (60 - jt') [ ( 60- jl^) “tOj^ ] 



1 - tOp^ (tt> - jt^) 

60 [ ( 60- jv)^-C^)^] * 

?= “b 



60[( 60- jp)2 - 60 j^2 ] 



, Z = 1 - 



to, 



60 ( 60 - ]V ) 



and substitute into (62) to obtain, 

VxH = a (EX-jE ?) + a (E X+jE Y) + a E Z. 

3 toco XX y yy x 



(62) 



(63) 

(64) 



(65) 
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For cross reference to literature employing the standard URSI symbols, 



X = Y - r = 1 - j u/ Oi . 



03' 



_k 

03 



the symbols in (63) and (64) are equivalent to, 

X = 1- X r , Y = X Y , Z = 1 - 2L 

r 



( 66 ) 



(67) 



2 2 

r - Y 



2 2 

r - Y 



The symbols (67) will be used later, but for the present it is more conven- 
ient to retain the symbols (63) and (64). 



3.32 Rotating Coordinates 

The a and a components of equation (65) will be broken into left 

X y 

and right hand circular components. Re-arranging equfi tion (65), 



1 



J 03C. 



vxH=a (E X-jE Y) + ja (EY - jE X) + a E Z, 



X X 



y X 



z z 



or. 



vxH= E (a X+ja Y)-jE (a Y +ja X) 



] 03Co A A -r A 

Now let, 

a X + ja Y= A, (a +ja )+A_(a - ja) , 

X y 1 x y 2'x y 

and solve for the undetermined constants A^ and A^ , 



Ai + Aj=X 

Ai-A2=Y 



r>U ^ 

X + Y 



A, = 

1 2 

A =1-:^ 
2 2 



also let. 



a Y+ja X — B, (a +ja ) + B_ (a “ j 3 ) 

X y lx y 2 X y 



( 68 ) 

(69) 



(70) 



(71) 



(72) 



and obtain, 



rsu 



Bi+B2=Y 



Bi-B2=X 



P - X +- Y = a 
®1 2 — 1 

B =X:i^=-A. 
2 2 2 



Substituting from (71) and (73) into equation (69) 

, X + Y / - 



J 03Co 



V X H = (E - j E ) 
X y 



(a + ja )+(E + jE)x 
X y X y 



X-Y(a -ia)+£EZ. 
— T~ X y 2 z 



(73) 



( 74 ) 
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From equations (63) and (64), 

X ±Y= 1- tOfa COp^ 

w[ ( to- ji;)^ - t 0 j^ 2 ] 

1 - X 

r ± Y ' 



O 

to [(to- jv) T to. ] 

1 •” — £1 ^ 

to[(to- iu) + tOj^][(to- iu) - tOj^] 

(75) 



upon letting, 

E = 1/2 Ve^ + E ^ . (h = tan"^ . (76) 

XY X y XY E 

X 



equation (74) may be written. 



7 X H = j 



w e 



E e 
XY 



XY 



( 1 - X ) ( a - j aM + a E_ ( 1 - ) 

P_ Y X y 2 z 



r + Y' • X 

p.) ! 



a + j a ) + E 

X y 






XY 



XY 



(77) 



3.321 The Complex Dielectric Tensor 

Designating the left hand rotating unit vector by a and the right hand 

J-i 

vector by a , that is , 

R 

a_ = a + j a , a„ = a - j a , (78) 

Lx y R X y 



and letting. 



€ / = 1 - X = 1 - X , e/=l-_2L. 

r+Y r-Y r 



V X H = j 00 Co 



(79) 



les , 


- j 0 1 




' ^ 


r c' 0 o~i 




E e 






L 




XY 




L 


o 

o 


• 


r. 3 0 


# 


- 


K 

0 0 Co' 




E e ^ XY 

XY 










E 




a 




z 




z 



= jto D. (80) 



3.322 The Complex Indices of Refraction 

Equation (80) yields three characteristic waves propagating in a 
homogeneous anisotropic ionosphere. There are two circularly polarized 
waves within the XY - plane rotating in opposite senses, and a linear 
wave parallel to the Z-axis. These waves have three distinct indices 
of refraction. 
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1#^ «4i» 



1 




I 



I 



( 81 ) 



Rationalizing the denominators of equation (79), 





IT 0) 2 n 

1 - P 0 


- J 


1 

CM 

a 

3 

1 


1 


co2 + p 
2 1 


■ 


_Ci)(co2 + j^2)^ 
\ 




CO 



P ^b) 



2 2 

C0[(C0 + CO]^) + V 1 

COp^ (CO - CO b) 

CO [(cO - COj^)^ + t'^] 



“ j 



V 






- J 



C0[(C0 + CO]^) 

2 

_CO [ (cO - COj^)' 




(82) 



(83) 



The corresponding indices of refraction will be obtained by extracting the 
square roots of equations (81), (82) and (83). The ordinary index of refrac- 
tion, n^ , is identical with that of equation (27). Upon applying equations 
(26) to equations (81) and (82), the left and right hand indices become. 



n 



L 

R 



= J_ 

n 



T . (m tOfe) 1 Jr _ gi (g, t 

,L co[(co ± CO, )^ + ] -J vL ± ^ 



)^+ P^]) 



2 2 



1 



1 / 2 ' 



CO 



[(co±o}^)2h- p2] 



■/ fl - “o' " 

\ [^ CO [(cO ± CO 



1/2 

7 T 



[(cO 

- I 1. CO 2(w±^ 



^^Ib 

^2 . 2. 

k) + 3 



2 2 
COp P 



CO r(co ± cot))^+ 3 

J 



2 2 

CO [ (co ±cOj^) + p 3 




1/2 

(84) 



It should be kept in mind that equation (80) is normalized, and hence 
indices (84) apply only for propagation of these characteristic waves. For 
wave combinations other than the characteristic waves, there will be 
coupling elements within the dielectric tensor, and a dispersion equation 
will be required for determining the indices. Which of these waves appear 
depends upon the direction of incidence^as well as upon the frequency. 



3.323 Physical Interpretation 

The critical angles of reflection, corresponding to the real parts of 
the indices of refraction in equation (84), may be determined by equation 
(30). This will be done for a low loss ionosphere, p « 0, and a physical 
interpretation formulated. 
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Setting = 0 in (84) and squaring, is equivalent to setting = 0 in 
equations (82) and (83). Hence, for applying equation (30), set 
2 2 



00 



00 ( CO ± OO]^) 



= sin 00 , 



2 2 
00 ± oOj^ oo - oOp^ sec 0 o= 0 . 



( 85 ) 



( 86 ) 



Taking the positive sign in equation ( 86 ) for oo^^^ , 



00 - = - + V ( 00 ^ sec 0 o)^ + 



cL 






(87) 



and taking the negative sign for o) 



cR ' 



oo^j^ = + + V (cOp sec Ql)^ + 

2 ’2 



( 88 ) 



Since co determined by equation (81) is identical with oo determined 

-w,. ^ 

by equation (33) , 



m = oo„ = oo„ sec 0 , 



while , 



I TT"^ ^ oo , 

"cL = ’'“o + ' 

‘^cR ■ ^ ^c ^ W 1 ^ 




2 2 

For frequencies sufficiently high such that (oii ) >> {<^h] 



^c 

Hence , 
00 



00 x2 



1 

l2 «) 









1+ 






8 0)^ 



— 60 . (90) 

^ 8 CO 

c 



cL 



+ 



V - "b 



8 oo. 



« 00, 



(89) 



(91) 



and. 



oo_r. « oo + « oi + ^b 



"cR 



8 00 2 
c 



(92) 



Since the three characteristic waves have distinct indices of refraction, 
they propagate within the ionosphere along different ray paths, and have 
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different attenuations. The attenuations spoken of above are those due 
to C n] . There are other attenuations due to numerous anamolies of 
the ionosphere. Because of the different critical frequencies, the right 
hand wave may be reflected at a frequency for which the other waves 
pass through the ionosphere. 

Recalling that the series in equation (54) diverges for frequencies 
less than the cyclotron frequency (or gyro frequency), it may be con- 




the extraordinary wave confronts a conducting boundary and is reflected 

as in path c (Fig. 6). For frequencies such that to, < to < co , the 
~ be 

wave enters the ionosphere, suffers a refraction, and is refracted back 
to the earth, as in path b (Fig. 6). On the other hand, for frequencies 
to > to D / the wave is slightly refracted but passes through the 
ionosphere as in path ^ (Fig. 6). The three waves do not necessarily 
all exist simultaneously, depending upon the angle of the ray path 
with the geodesic field. 

3.4 Azimuth Component of Exciting Field 

In the preceding analysis, the exciting field was considered 
vertically polarized. However, a wave incident upon the ionosphere 
from a horizontally polarized antenna will have a horizontal component 
in addition to a vertically polarized component. Hence, it becomes 
necessary to consider the azimuth component of an exciting field. 

Re-writing equation (3 6) for convenience, 

F = -e E = m d V + pmV+eVxB , t = -NeV , (36) 

dt ^ 

and considering the exciting field in the form. 
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(93) 



= - eVo X B= j m (o3 - iv) 



E = E„ (- a sin(p+a cos <o) , 

<P X z 

the iteration follows analogously to that for E . Accordingly, 

0 

following steps (43), etc., 

Vo = - e Q (0 
j m (o3- jv) 

V , = e2 E„ a C5 ^ ® ^o 05^ (-a sin cp+ a^ cos d x a^ ^ 

[J m(o}- Jp) ] m (co - jp)^ 

e Eo OJj^ sin <p a 
£ ! X — , 

m (co - jv)^ 

7 ~ ~ 

- e E„ B cOu sin (p a x a . „ 2 • 

V 2 = - !z b ^ X z_ = -leE o sin <p a 

j { u)- j vr — siv-'jV)3 — 



- ^ j e2 E„B^ oJb^ sin (pa^x a^ ^ eE^ t0j^3 sin <p a 
^ j m2 ( CO - iv 

—2 E, B^ cOb^ ^ 

j m 2 ( ^ - ji^)5 m (oJ - ]v)^ 



^4 = 



-b jv , 

m ( to - ivr 

■e2 E„ B^ cOb^ sin (pa^x a^ ^ -j e E„ cOb"^ sin <pa^ 



Substituting the above into curl H , using l = -NeV , and 

2 



considering Eg a = a E _ + a_ E 

3 W €0 



o ^ ^ X, . ^ X, , curl H = 

<p y y z z 

E„a r 1- ] - ja_ E 



(P 



3 5 

to, av 



to (a' - jp) 



X y 



to [(to - 



ii^)' 



_^b + ^b + - 

(a? - ]v)^ (a' - iv )^ 



av 



-]-v,v[ 

1 } 






(a' - 31 ^)^ 



(94) 



(95) 

(96) 

(97) 

(98) 



(99) 



(w - jp)^ TaT-jiTp 
Summing the series over the range of frequencies for which the series 
converges , 

curl H = j a' €c 






1- ~ | -ja E % % 

oi a'- iv) J ^ ^ a (a* - iv)^ 



1 


2 2 
-i E ‘"b 


1 — • 


1_ w b''^ 




(a - jp)^ 


^ ^ a (a - jp)^ 


1 

in: 

1 



(100) 
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or. 



curl H= j toCo ^ 

. 

[a - j a ^'b 

X y — 



E, [ 1- 

^ . a' (o!- ji^) 



] - jE ^ 

^ 6D[ (to - jp)2 - OJ, ^ ] 



] 

ui - iv) ^ 

Substituting for a E. of equation (93), and using E and E , 



(101) 



9 

curl H = j a'€. 



E O) 6lV 
“J ^ V P b 



2 2 
CO ^ CO, 
p b 



C0[( 05 - jl^)^ - CO, 



+ a E J 1- ^P 

y y 



a' (co- ju) 



+aE £1-^P 
z z ^ 



CO (05 - jp) t( CO- jp) " COj^^^ 



60 ( 05 - jv)‘ 



(102) 



Let 



Substituting from equations (63) and (64), 
curl H = j 60 €„/-jE (a Y+ja X) + a E Z 

] y X y z z 
a Y+ja X=A, (a + ja ) + A_ (a — j 9 ) / 

X y 1 X y 2 X y 



(103) 

(104) 



and determine the arbitrary coefficients A^ and A^ , 

■ - 

Aj-Aj.X . 

Thus, (104) becomes. 



A^ = 1/2 (Y + X) 
A^ = 1/2 (Y - X) . 



(105) 



^ curl H = - j E ( X + y) (a + j a )+ j E (X - y) X 
J y X y y 

2 2 

(a - j a ) + a E z . 

X y z z 



(106) 



Changing to the URSI symbols by substituting from equations (67), (75) 
and (78) , 

curl H= j o5€o J - J \ (1- X ) a + j E (1- X ) a + 

-J r + Y ^ r-Y 



E (1- XJ a 

" r ^ 

Using equations (79), _ 



Curl H = j 05 €, 



'L 

0 € 



(107) 



0 0 

R 0 

0 e: 



-> E/2 

J E/2 



R 



= j CO D 



a 08) 
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Hence, the same characteristic waves appear as in the case of the 
vertically polarized exciting wave. Which of these waves appear 
depends upon the angle of incidence and direction of incidence, 
as well as upon the frequency. 



3.5 Horizontally polarized exciting field 

Suppose an exciting electric field is incident upon the ionosphere 
from a horizontally polarized antenna, the incident field being in the 
form , 

j a' t - V . r 



(109) 



- a^ being shown in figure 4. It may be seen that - a^ is composed of 

two orthogonal components in the a and a directions. 

a <p 

If, in equations (40), is replaced with , and if in equation 
(93), E„ is replaced with E- , then from table (35), E.of equation (109) 



may be expressed by. 









with. 



E = E COS_0_jin_^^^ 

^ 1 - sin2 0 sin^ (p ' V 1- sin^ 0 sin^ ^ 



Also, in equations (40), 



/~p 2 , p 2 = Eo r 2 . . 2 

V E + b V QQ3 0sin <p , 



and, 



tan 0 



XV 



= y = - cot 0 cos (p . 



(110) 

( 111 ) 

(112) 

(113) 



Therefore, for the horizontally polarized exciting electric field vector, 
equations (77) and (107) may be combined intdt 

{> 

+ a^ Co cos 0 sin <p cos 0 sin <p E^, 



curl H = j coCo -I 1/2 [ e,. ' \y~ aj^+ e XV Cj^ ] J ^ 



J \ - sin^ 



,, . 2 
0s in (p 



sin 0COS ^ 
(114) 






+ jWfoj J 1/2 ^ ^z 1 °^ ^ ? 

L _J y 1" sin'^ 0 sin^<p 

To determine the critical frequencies involved in equation (114), it would 
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seem preferable to consider individually the various normalized components. 



4. THE DISPERSION EQUATION 

4.1 Maxwell's Equations in Complex Form 

The instantaneous iy;fexwell equations, 

div B = 0, div D= p , curl E = - 5B , curl H = d D 

9t St' 

L=aE=pV = - NeV ,D = €'CoE ,B = /Lt'poH , 



are possibly more useful for steady state time harmonic cases when con- 
sidered in the form, 

div B = Ojdiv D = 0, curl E = - j w jliH, curl H = t+ j W C E . (116) 

To show that div D = 0, first formulate the equation for the continuity 
of charge by taking div curl H , 



div curl H = 0 = 


div 


L + 9 div D = div l + 9 p 
9t 9t ' 


(117) 


or 








div L = - 9p 

TT . 






(118) 


But, for conductors. 








div i = crdiv E = 


cr 


div D = cr p , 


(119) 


££-+ 2_ p= 0 


e 


€ 


(120) 



9t c 



The solution of (120) may be written 

p^p^e'lf/oU . (121) 

From equation (121), a charge density within a lossy medium must vanish 
with time, and hence must vanish within a steady state condition. Of 
course, since a dielectric is considered to be free of charges, div D also 
vanishes within a dielectric. 

Now consider the case of a spherical wave front given by. 



E = Ep(e, (p) e. 



j 03t - y . r 



( 122 ) 



which is sufficiently remote from the source for. 



9 E, 



0 , 9 Eo « 0 , 



to hold locally. Also assume y = a^ y . Then, 

div Q~y * ^ ) = 1 r_^ 

sin 9/ 9r L ^ 



Eo r sin e y . r - 



-a F j " f • ^ , ~ y • ^ I "“E 

— a^ . Lq ^ — yr Q ' +p> ' t= — v.'*^ 



y • ° e 

r 



. a 



- y . r 



} 



i- r . i, e 

r2 



(123) 



But, the wave is assumed to be a radiated wave such that the inverse 
square terms vanish. Hence, in practice, for a radiated wave. 



div 



Eo (6 , (p ) ^ - y . r 



- “ y • En (6,1 o) Q 

r 



- y . r 



Likewise, consider 
curl (_Ej 3 _ e ^ ‘ ^) = -li 



sin e 



r 

ydr 



r a 



6 

d/B0 



IjilfT y- ^ rE 



=0 e- r • 



(124) 



r sin 0 a 

B/d(p 



(P 



- y . r 



LSln^j; e 

COS 0 - y. r ^ 

= a sin 0 ® - a r sin 0 , 4 • r ^ ^ 

r e Eo (-re ) + 

r sin 0 ^ 






) . 



(125) 



Again, assuming a radiated field for which the inverse square termsvanish. 



curl ( E^e “ ‘ = a^. 0 - a^ ( - y E 



- y • r ) • 



0 <p e* * '') + " r Eo e 

r r 



(126) 



But, 



y^ Eo = 



y 

E„ 



<p 



0 g Eo 



<P 



= a 0-a.yEo +a y Eo_ 
r 9 0 (f> e 



(127) 
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Therefore, in practice, 
curl (8 . (d) g 

Also, in case the relative constant is a tensor, e', then 
D = . E . 



-y 






X Eq (6 , <p) 



- y 



(128) 

(129) 



Therefore, the complex form of Maxwell's equations for radiated fields 
may be written, 

y . H = 0, y . €oC'. E = ,0, l = - NeV, y x 2 i= j co/Jo H, - y x H = 

-NeV + j W€oE = j w €o€ ' . E , (130) 

in which, e ^is defined by 

"e ' . E = - NeV/ jwCo + E , (131) 

The wave equation becomes, 

y X ( y X E) - k?o e' . E = 0 . (132) 



4. 11 Orientation of field Vectors 

From equations (129) and (13 0), the relative orientation of the field 
vectors may be obtained. For this purpose, consider 

y . 5 = 0,"y . H = 0 , (133) 

from which it may be concluded that y is normal to both D and H. 

Then consider, 

-yxH = +jcoD, (134) 

from which it appears that D is normal to bothy and H. Therefore, the 
three vectors D, H, and y constitute an orthogonal triple. Finally 
consider 

y X E = j 60 Mo ii / (135) 

from which it is seen that H is also normal to the plane containing yand E. 

Assuming y to bq^t an arbitrary angle ij) with the geodesic field Bo , 
the relative orientation of the unit vectors is then illustrated in figure 7. 
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Thus, the TEM incident wave becomes a complex TM wave within the 
ionosphere. 



4.2 The Dispersion Equation 

For y = j ko n , the wave equation reduces to, 

n X (n X E ) + e ' . E = 0 . (13 6) 

upon expanding the cross products and providing for all components of the 
dielectric tensor, equation (133) may be expressed as. 



, 2 2 , 

-(n + n ) 

y z 



n n 

y X 

n n 
z X 



n n 

X y 



n n 
z y 



= 0 



Combining the tensors , 



n n 
X z 



-(n^ + n n n 
V z v z 



^ 2 ^ 2 . 

-(n + n ) 

X y 





"e 




P €' e' 




X 




XX xy xa. 


• 


E 

y 


+ 


c' c' c' 

y X yy yz 




E 








z 




c zy zz 








z X 








- 



€ ' 


-n 2 -n 2 


( + 


n n 


/ 

e 


+ 


n 


n 






XX 


y z 


xy 


X 


y 


Xz 




X 


z 






e' 


+n n 


f ' - 


n2 -n2 


c' 


+ 


n 


n 






y 


y X 


yy 


X 


Z 


yz 




y 


z 


• 




e' 


+n n 


e' + 


n 


n 


c' 


_ 


n 2 


2 

•^n 






z X 


z X 


zy 


2 


y 


zz 




y 


X 







= 0 



(13 7) 



(138) 



Equations (136) constitute three homogeneous linear equations defin- 
ing the components of E. A necessary and sufficient condition that a 
solution other than the trivial exists, is that the determinant of the 
coefficient matrix vanishes.. The resulting characteristic equation is 
called the dispersion equation, as its eigenvalues determine the indices 
of refraction and coefficients of attenuation. That is, the dispersion 
equation is , 



€' 


•n 2 - n 2 


e ' 


X X 


y z 


xy 


c' + 


n n 


i 


y X 


y X 


'yy 


e' + 


n n 




Z X 


z X 


zy 



+ n n 
X y 



€ + n n 

XZ X z 



2 2 e ' L 

-n^-n yz+ y z 
X z 

+ n n 

z y 



= 0 



(139) 



, 2 2 
€ - n — n 

zz y X 
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4.21 Special Cases 

As an indication of the application of the dispersion equation, some 

special cases will be considered. For this purpose, re-write equations 

(68) and (103), respectively, in the following forms. For E , 

0 



- y X H = j co€ 



1- X r -j X Y 

p2 _ - Y 2 

j X Y 1 - X r 

_ y 2 p2 _ y2 



0 

0 

1 - 

r 



and for E 



<P 

- y X H = j U) Co 



Let 



(140) 



0 -j X Y 

-y2 

0 -1- X r 
0 0 



0 

1- X 



(141) 



y = ya = j k^n = j kon (a cos 0+ a sin 0 cos <p+ a sin 0 cos (p) (142) 

X Y ^ 



Using the dielectric tensor from equation (140) along with the dispersion 
equation (139), and also substituting from equation (142), for (b = tt/2. 



1- r? sin^ 6 - TX - j X Y 

pi ^ 

^ 1 - n^ - rx ^ 

_ y2 p2 - Y^ 

2 

n sin 0COS 0 0 



n sin 0 cos 0 



1-n^ cos^ 0 - X 



= 0 (143) 



and for <p = 0, 

1 “ rx _ „2 2 0 ~ 1 ^ — + n^ sin 0 cos 0 



r2 - Y^ 



p2 _ x2 



1 X Y ^ 2 . ^ ^ 

^72 ^2 0COS 0 1- rx 

^ ■ r2 - y2 



- n2 



n^ cos ^ 0 



1- X 



-n 



=0 



(144) 
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Equation (143) is the dispersion equation for a vertically polarized 
wave propagating in a longitudinal direction, whereas equation (144) 
is for a vertically polarized wave propagating transverse to the geodesic 



field . 



Now, in equation (143), let 6 = tt/2, 

1- r? - T X -i X Y 

p2 _ y2 r ^ 



i X Y 



2 2 
- x^ 



l-n -rj(. 



0 



r2 _ x2 

0 1- r X 



- x^ 



= 0 



Thus , 



( 1 - r X - n^)^ - fx Y )^ = 0, 



- Y^ 



" Y" 



and 



n = 1 - X 



r ± Y 

For equation (144) with 0 = tt/2, 

1 - - r X -i X Y 

r ^ - Y^ - Y^ 



\ X Y 
r'-Y^ 



1 - rx 



r2 - y2 
0 



1- X 



- n 



= 0 . 



one solution is , 
n^ = 1 - X 



For the other solution. 



(1 - X -n^) (1- r X ) ‘ ( X Y ) - 0, 

„“2 Z 

r - Y 



- Y^ 



2 2 

r - Y 



n = 1 - r X 

- y2 



2 2 
X Y 



(r^ - rx - Y^) (r^ - y^) , 



or, 2 

n = 1 - X_ 

r - Y 



r- X 



(145) 

(146) 

(147) 

(148) 

(149) 

(150) 

(151) 

(152) 
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Solution (147) is for a vertically polarized wave in the longitudinal 
direction whose exciting field is normal to the geodesic field, and 
solutions (149) and (152) are analogous solutions for the transverse case. 

For the azimuth, or horizontal component of the exciting field, the 
dispersion equation for <p = IT /2, that is, for longitudinal propagation, 
becomes , 



2 ■ 2 ,, 
-n^ sin 6 



-j XY 

r'-Y^ 
1-n^ - r X 



2 

n sin 0 cos 0 



- X^ 



n sin 0 cos 0 



l-n^ cos 2 0 - X 

r 



= 0 



(153) 



Expanding, 



i X Y (0) + (1- r X 



2 2 

r - Y 



- x^ 



- n 



2 ) 



from which a solution is obtained, 
n^ = 1 - r X 

j^2 _ y2 

For the other factor, 



- sin 0 n sin 0 cos 0 
cos 0 l-n2 cos 2 0- X 



= 0, (154) 



(155) 



-(1- X ) sin 0 + n sin 0 cos 0 - n sin 0 cos 0= -(1- X sin 0), 

r r 

which is constant. Therefore, solution (155) is the sole solution for this 



case. 



For transverse propagation, <p= 0, and the exciting field is parallel 
with the geodesic field. Hence, in this case, the ordinary solution 



holds, namely 



n^ = 1- 



X 



(156) 



5. APPLETON EQUATION FOR ARBITRARILY ORIENTED GEODESIC FIELD 
5.1 The Appleton Euqation 

The Appleton equation, sometimes referred to as the Appleton-Hartree 
equation is an equation for determining the complex index of refraction in 
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a homogeneous, anisotropic ionosphere. It is customarily derived by 
choosing one of tWo coordinate planes as being determined by the y 
and Bo vectors. It will be derived herein with an arbitrarily oriented 
geodesic field so that a greater leeway in the sfelection of the 
coordinate axes is permissible for applications. Refer to figure 8. 





The following symbols will be used , 

60 , = e . 60 ^ = Ne^ ,r = l“ii^/60,X= 60 ^ ,Y = ^b_ - 
m Pm c, ^0, 

X y z 



(157) 



Y = Y cos 0 , Y = Y sin 0 cos <p , Y = Y sin 0 sin (O = 
X y ^ z 



Y cos 0 = Yj^ , 


(158) 


Y^ = Y sin = Y 1- sin^ 0 sin ^ <p , e Bo = - jy . 

j 60 m 


(159) 


The geodesic field is taken as. 




Bo = Bo (a^cos 6+ sin 0 cos cp + a^ sin 0 sin tp ) , 


(160) 


and y as , 




y = a y = j k„ n . 


(161) 



5.2 Lorentz Conductivity Tensor for An Exciting Wave 

The exciting field will be postulated to vary as J e X • ^ • 

r 

Hence, Maxwell's equations may be written in the form, 

y . H = 0, y . foe E = 0, y X E = jco/ioH , - yx H = l + 
j 60 Co E, I = - N,eV . (162) 
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Solving for t » 

L = - y . X H - jo5€o E = - [ j coCo E + _1 yx(yxE) ] = - ja'c o A 

j W/io 

(163y 



[ E + - k ^ n X (n X E ) ] 



- k 



Now, 



or. 



nx(nxE) = nxfhx(a E +a E +a E)] = n'^a x(a E - 
^ X X Y y z z z y X 

a E ) , 

X y 

nx(nxE) = -n^(a E +a E +a 0). 

X X y y z 



(164) ' 

(165) 



Therefore , 



t = - j C 0 €o [ E + n X ( n X E)?= - j wCo C E - (a E + a E + a 0) ] , 

J X X y y z ' ^ ' 

(166) 

and hence, 

L = j 0)C 0 

(167) 



'n^-l 


0 


0 




1: 


0 


n2-l 


0 




X 


0 


0 


-1 




E 

y 






E 










Z ___ 



The Lorentz conductivity, a , is defined by, 

I = a . E . (168) 

2 2 2 

Hence, if M is defined by M = n - 1, the Lorentz conductivity tensor 
becomes , 

^ = =jo)„ 'A k;2 



n^ -1 Q 


o" 






O 


O 


O n -1 


O 




O 


m2 


O 


O O 


-1 


* 


O 


O 


-1 



M 



1 O O 

O 1 O 

0 0-1 



M^J 



(169) 



A corresponding Lorentz resistivity tensor, p , may be defined by, 
p . (, = p . O' . E = I . E = E , (170) 

witij I being the identity tenser. 
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The resistivity tensor may be found from equation (36), which will be 
repeated here for convenience, 

m d V + u m V+eV xBg= - eE. (36) 

d t 



In the steady state phasor form, equation (36) is, 
-ei = (jojm + i/m) V + eVxB, , 

Using equation (160) for B in V x B, , 



VxB, = 



or, 



a a a 

■ X y 2 

’ V V V 

X y 2 

cos 0 sin 6 cos <p sin 6sin^ 



V X B„= B.[ a (V sin0sin(£)-V sin0cos(p) + a (V cos 6 - V 
°® y y z 'y z ) 

sin0sin<p) + a (V sin0 cos <p - V cos 0 ) ] . 

2 X y 



(171) 



(172) 



(173) 



Substituting equation (173) into equation (171), factoring j o) m from 
the first two terms of (171) and using the symbol F, the resulting equation 
is, 

-e E = j turn r(a V + a V.*+a V ) + eB.[a (V sinSsincp- 
xxyy2z °xy 

V sin 0 cos <p ) + a (V cos 0-V sin0sin<p) + a (V sin 0 cos 0 

2 y z X z X 

- V cos 0 ) ] . 

y 

Substituting symbols (157) into the re-arrangement of equation (174) as 
a tensor equation, the equation may be written. 



= - 1 03 m 
e 



r 



-JY 



jY 



z - y 
r -jY 



-jY. jY, 





V 




X 


• 


V 




y 




V 




Z { 



= -1 



j 03 Cc,X 



r 

^^z 

-jY 



r 

jY 



jY 



-jY. 



(175) 



with the velocity matrix converted to the current density matrix by 
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multiplying and dividing by - Ne. 

From equations (170) and (175), the resistivity tensor may be taken 
to be , _ 



P 



- _1 

j OiSo ^ 



r 








T 








T 



( 176 ) 



Consider equation (170), and into it, substitute from equation (108). 

The result is , 

E = p.t = p.CT. E,I.E“p.a .E=0, (177) 

or, = = = 

[p.(7-l].E=0 . (178) 

Equations (178) constitute a set of the linearly homogeneous equations 
for the three components of E. Hence, for a non-trivial solution to exist, 

p . a -1 =0 . (179) 

Multiplying (169) by (176) and substituting into equation (179), the 

characteristic equation becomes, 

- 

X3 

m2 

(180) 



or. 



(r+ X )^ 

M^ 

( r+ 

^ m2 

(r+ X 2 Y ^ 
M^ > 



r+ 


- jY 


JY 




m2 


z 


y_ 

m2 






r + 


-JY 


= 0 


z 




_Y 

m2 




-jY 


JY 


r-x 




y X 

sterm inant. 


M^ 




' (r-x) 


+ jY Y 


Y - jY 


Y Y 


X y 


z X 


y : 


. ) -y| 

1 


( r - X) = 


^ 0, 





M 



r - X 



( r + i^) - = 0, 

M ^ 



(181) 

(182) 



2 2 2 

in which substitutions = Y and Y„ = Y„. + Y were made. 

L z T y 
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I 



Solving for F+ X 
M 

r + 2^ =1/2 



± J Y„ 



r -X 



( r - x2) 



+ 4 Y, 



J 



or, transposing. 



M 



Y 2 
r- T 



Y^ 

T 



2 (r- X) 



4 ( r- 



+ 4 Y, 



Taking the reciprocal, multiplying by X, and eliminating M , 

X 



i./t 2 2 , 

M = n -1 = - 



r- T 



Y„ 



+ Y, 



l4-l^)2 



(18 3) 



(184) 



(185) 



Solving for n and substituting from equations (157) and (159), 



2 

n 



= 1 - 





COS ?/) 



(0p2 / 

2(1-0) 2/ - jy/ to) 

1i72 



± 



o), ^/o)^ sin^j/) 
b + 

4(1- o) Vo)^ - jv Top 

p 

(186) 



Equation (186) is the well known Appleton equation with 0 being the angle 
between the geodesic field and the direction of propagation. It is not 
restricted to any coordinate system. That is, it is a mathematical 
invariant. 



6 . FARADAY ROTATION 

6.1 Waves Through the Ionosphere 

Waves passing through the ionosphere are, generally speaking, broken 
into two or three distinct waves having distinct indices of refraction, and 
they travel by distinct ray paths. If the waves are attenuated, they are 
unlikely to recombine into linearly polarized waves. 
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However, if attenuation is negligible, the circularly polarized 
waves may recombine into linearly polarized waves, but they will have 
their plane of polarization rotated from the original plane. This is 
referred to as Faraday rotation. 

To examine this rotation analytically, consider two such waves 
emerging from the ionosphere after undergoing different path length 



shifts in phases, (t)^, <t>^ , respectively. Write the phasor equations, 
-Jkofo -jk„ .fn , • ds , ^ -jk, Jn 2 . di^ J 



E - Eo e 



Since the ionosphere is homogeneous by hypothesis, 
E = E. e (e "l "l + e ’i’'” "2 ^2 ) 

or, ,_ . n/s, + n,, s 



E = E, 



11 



!i_ 2 ) 



<"l=l - " 2 " 2 > ] 



[• 



e^<Vl-"2 =2> + 



(187) 

(188) 



(189) 



in which s^ and s^ are the respective path lengths. 

Upon multiplying and dividing by 2 , 

E = 2 Eo cos [ ^ (n s - n s J ] ^^ 1^1 + ^2®2^ 

2 ^ ^ ^ ^ 

Thus, the resultant of the shifted vector is 2EoCOS [.ko(ns - ns ], 

2 ^ ^ 2 2 

J<g_(n^Sj - ^ 1282 ) is the angle of the resultant. This may be verified 
by referring to figure 9 . 

<Pl <p2 

? 

<Pi = ^0 ns 

^2 ^ ^ ”2 ®2 
2 



From equation (190), it may be seen that the phase shift of the field 
through the medium is determined by the average path length of the two 
ray paths. 
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i 



c 



1 




« 

A 



i 






6 . 2 Polarization 



In order to determine a measure of the complex polarization of a 
wave propagating within a homogeneous anisotropic ionosphere, it is 
desirable to find the ratios of the electric field components as deter- 
mined by equation (179). The matrix form is. 



M 



X 



r + X 
« — 




- jY, 


M^ 






jY. 


r+ 


X 


z 






-jY 

y 







-JY 



M 



JY 



X_ 
m2 



- (r- X) 






L 



= 0 . 



From the theory of linear homogeneous equations, the ratio E : 

X 

Ey : E^ may be found from either pair of the three equations, given by 
equation (191), by omitting the first, second, and third columns of 
coefficients, respectively. The signs are alternately plus , minus, 
plus . Thus , 



E = 
z 


r + X 

m2 


jY 




jY jY 

m2 




jY 


-(r - X) 




1 

1 




X 9 

M^ 




Y 9 

m2 



-JY 



r+ X 

M 
jY.. 



The polarization R is usually defined in terms of the ratio of two 
components of the electric field normal to the direction of --phase pro- 
pagation. To facilitate the algebra, let, 

Y = 0 , Y = , and Y = Y^ 

X z L y T 

This, in effect, rotates the coordinate axes (figure 8) such that B, 
lies in the Yz - plane. Now write the ratio, 

R = E 



r+ X jY 
"72 X 






jY. 



“(r - X) 



(r-x)(r+x )-Y 



M' 



-JY 



-(r- X) 



] (r - X) Y + Y Y 
z X y 



(191) 



(192) 



(193) 



(194) 
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which, upon substituting from equations (183) and (193), becomes 

“3 ; 



R = _i 



f"Y 2 



2Y. 



X 



i(r- x)^ 



+ 4Y 



= -1 

JY, 



(r + j^) 



The other ratios may be expressed. 




-jY^ 

r - X 



E 

z 



E 

y 



R. 

r - X 



( 195 ) 



(196) 



7. CONCLUSION 

7 . 1 Equivalence of Two Points of View 

For the purpose of yielding a better insight into the mechanics involved 
when a wave is propagated into the ionosphere, the principles of electron 
ballistics were applied to the free electrons. In particular, equation (36) 
was solved by an iterative procedure^ in the procedure . for deriving 
equation (68), the latter equation yielding the dielectric tensor. 

The components of the dielectric tensor were used in equation (139) 
for determining the indices of refraction. This dispersion equation can be 
formally written, 

[n^r-(n +“€')].§ = 0, (197) 

in which an index tensor n is introduced, with the components of n 
being defined as , 

n = [ n.. ] = [n.i n. ] ,i, j = x»y, z. (198) 

Equation (36) was also used as a key equation of constraint in 
deriving the Appleton equation (186). Thus equation (36) may be thought 
of as a sort of common denominator between the two procedures for 
finding the complex indices of refraction. 

In fact, equation (68) can be derived much more compactly by 
formalized procedures. For this purpose, consider the ccj’re spending 
complex form of equation (36), 

-eE=(jo)ni+pm)V + eVxB^ (199) 



36 



Referring to figure 5, and substituting 



-JY = e B, 

— / 

j 05 m 

equation (199) becomes, 

-g E = j 05 m [( 1 - j V /oi)y - \ Y V X a^ ] = j 05 m x 

[rV+jYVa -jYVa] 

X y y X 



( 200 ) 



(2 01 ) 



Hence, 
E 



X 



E 

L z 



05 m 
e 



■]Y 0 



jY r 

0 0 



0 

r 



V 



V 



V 



= -1 



Co X 



r -jY 0 

jY r 0 

0 or 



( 202 ) 



From equation (202), the Lorentz resistivity tensor p may be written 
namely, 

' r -jY 0 

jY r 0 

0 0 r 



= - 1 
j05 CoX 



(203) 



Now, the complex dielectric tensor is given by, 
= ' = I + j o5€o (p)-l 



(204) 



Thus, the inverse resistivity tensor may be found by customary matrix 



algebra , 

^ -1 = “ja>€ X 



r(r^ - X^ ) 



or , = , 

P -1 = -] 05 Co 



r X 



r 

-jr Y r 

0 0 
iXY 



jPY 0 
2 






“iXY rx 

x2- ^ x^ -r^ 



0 

- Y^ J , 



0 

X 



(205) 



(206) 



3S 



Substituting into equation (204) and introducing symbols (63) and 



( 64 ), 

€ 















1- r X 


-jXY 






X 


-jY 


0 


x2 ^ p2 


x2 - 


0 






iXY 


1 - r X 


0 




JY 


X 


0 


x2 . p2 


x2 _ r 2 




















0 


0 


Z 


0 0 


0 


1- X 














r. 




— 




— 



( 207 ) 



Substituting e' into 

curl H = j o)€o e 'E , (208) 

yields precisely equation (68), that is, _1 curl H = 

j 

a (E X-JE Y) +a ( jE Y + E X) + a z" . (68) 

X X y y X y z 

The above derivation tends to place more confidence in the previous 
physical interpretations. It also serves to tie the iterative procedure to 
the Appleton equation. 



7 . 2 Inhomogeneous Anisotropic Ionosphere 

There is no exact mathematical model for the inhomogeneous aniso- 
tropic ionosphere. Many statistical measurements and mathematical 
interpretations have been and are being made. A vast amount of lit- 
erature exists, but various precise studies remain to be made. 

A detailed study of the inhomogeneous ionosphere was entirely 
beyond the scope of the time and fadilities available for the pre- 
paration of this report. 



UNCLASSIFIED 

Security Classification 



DOCUMENT CONTROL DATA • RAD 


(SmmifUr ctm90Hicmtion oi tittm, kody oi mbmirmct indmKind m\notmtion mvt bm mttmfd wdimn owrmit fpmrl i« 


1. ORIGINATING ACTIVITY (Corpormtm muthcr) 




U. S. Naval Postgraduate School 


UNCLASSIFIED 


Zb GROUP 


Monterey, California 


J. REPORT TITLE 




Electron Ballistics and Electromagnetic Waves in the Ionosphere 


4. DESCRIPTIVE NOTES (Typm oi rmport mnd tnchfirm dmtmm) 




Technical Report, 1965. 





Chaney, Jesse G. 



REPOHT DATE 

10 November 1965 



?«• TOTAL NO. OF 

44 



7 *. 



la. CONTRACT OR GRANT NO. 
b. PROJECT NO. 



Work Request No, 
62001 



f «. ORIGINATOR'S REPORT NUMBER(5; 

TR No. 56 



$b. OTHER NBRORT NOfS; (Any oUft numb»f tint mmy bu 
tftia nport) 



10. A VA IL ABILITY/LIMITATION NOTICES 

Unlimited 



I. SUPPLEMENTARY NOTES 



12. SPONSORINO MILITARY ACTIVITY 

United States Naval Missile Center 



IS. ABSTRACT 



It is shown that ion collisions introduce a slight reduction in the 
plasma frequency, along v\i th an exponential decay of transient electron 
oscillations. The critical frequencies for penetration of a homogeneous 
ionosphere, for both isotropic and anisotropic ionosphere^ are determined. 
The characteristic waves, for electromagnetic propagation within a homo- 
geneous anisotropic ionosphere, are developed by considering an infinite 
series of electron velocities, produced by an exciting electric field. The 
complex indices of refraction are determined, both from a dispersion 
equation and from a derivation of the Appleton equation, which uses an 
arbitrary selection of the coordinate axes, thus emphasizing the invariance 
of the Appleton equation. Vector and tensor algebra is used throughout the 
analytical developments. 



DD 



FORM 

1 JAN 64 



1473 



38 



■IIMOLASSiriED 



Security Classification 



UNCIASSIFIED 

Security Classification 



KEY WORDS 


LINK A 


LINK B 


LINK C 


ROLE 


WT 


ROLE 


WT 


ROLE 


WT 


Ionosphere 
Electron ballistics 
Electromagnetic waves 















INSTRUCTIONS 



1. ORIGINATING ACTIVITY: Enter the name and address 
of the contrecto/;, subcontractor, grantee, Department of De- 
fense activity or other organization ('corporate author) issuing 
the report. 

2*. REPORT SECURITY CLASSIFICATION: Enter the over- 
tU security classification of the report* Indicate whether 
’’Restricted Data'* is included. Marking is to be in accord- 
ance with appropriate security regulations. 

2b. GROUP; Automatic downgrading is specified in DoD Di- 
rective 5200. 10 and Armed Forces Industrial Manual. Enter 
the group nximber. Also, when applicable, show that optional 
markings have been used for Group 3 and Group 4 as author- 
ized. 

3. REPORT TITLE: Enter the complete report title in all 
capital letters. Titles in all cases should be unclassified. 

If a meaningful title cannot be selected without classifica- 
tion, show title classification in all capitals in parenthesis 
immediately following the title. 

4. DESCRIPTIVE NOTES: If appropriate, enter the type of 
raport, e.g., interim, progress, summary, annual, or fitial. 

Give the inclusive dates when a specific reporting period is 
covered. 

5. AUTHOR(S): Enter the name(s) of author(a) as shown on 
or in the report. Enter last name, first name, middle initial. 

If military, show rank and branch of service. The name of 
the principal author is an absolute minimum requirement. 

6. REPORT DATE: Enter the date of the report as day, 
month, year; or month, year. If more than one date appears 
on the report, uae date of publication. 

U, TOTAL NUMBER OF PAGES: The total page count 
should follow normal pagination procedures, i.e., enter the 
number of pages containing information. 

7b. NUMBER OF REFERENCES: Enter the total number of 
references cited in the report* 

8*. CONTRACT OR GRANT NUMBER: If appropriate, enter 
the applicable number of the contract or grant under which 
the report was written. 

8b, 8c, & 8cf. PROJECT NUMBER: Enter the appropriate 
military department identification, such as project number, 
subproject number, system numbers, task number, etc. 

9e. ORIGINATOR'S REPORT NUMBER(S): Enter the offi- 
cial report number by which the document will be identified 
snd controlled by the originating activity. 'This number must 
be unique to this report. 

9b. OTHER REPORT NUMBER(S): If the report has been 
-ttltgnad any other report numbers by the atiginwtor 

or by the sponsor), also enter this number(s). 

10. AVAILABILITY/LIMITATION NOTICES: Enter any lim- 
itations on further dissemination of the report, other than those 



imposed by security classification, using standard statements 
such as: 

(1) "Qualified requesters may obtain copies of this 
report from DDC 

(2) "Foreign announcement and dissemination of this 
report by DDC is not authorized. " 

(3) "U. S. Government agencies may obtain copies of 
this report directly from DDC. Other qualified DDC 
users shall request through 

tf 



(4) "U. S. military agencies may obtain copies of this 

report directly from DDC Other qualified uaera 
shall request through 



(5) "All distribution of this report is controlled. Qual- 
ified DDC usera shall request through 



If the report has been furnished to the Office of Technical 
Services, Department of Commerce, for sale to the ptiblic, indi- 
cate this fact and enter the price, if known. 

XL SUPPLEMENTARY NOTES: Uae for addiUonal explana- 
tory notes. 

12. SPONSORING MILITARY ACTIVITY: Enter the name of 
the departmental project office or laboratory sponsoring (psjr- 
ing for) the research and development. Include address. 

13. ABSTRACT: Enter an abstract giving a brief and factual 
summary of the document indicative of the report, even though 
it may also appear elsewhere in the body of the technical re- 
port. If additional space is required, a continuation sheet shall 
be attached. 

It is highly desirable that the abstract of claaaified reports 
be unclassified. Each paragraph of the abstract shall end with 
an indication of the military security classification of the in- 
formation in the paragraph, represented as (TS), (S), (C), or (U). 

There is no limitation on the length of the abstract. How- 
ever, the suggested length is from 150 to 225 words. 

14. KEY WORDS: Key words are technically meaningful terms 
or short phrases that characterize a report and may be used as 
index entries for cataloging the report. Key words must be 
selected so that no security clasaification is required. Identi- 
fiers, such as equipment model designation, trade nasM, military 
project code name, geographic location, siay be used as key 
words but will be followed by an indication of technical con- 
tut. The assignment of Unl^, r«Us, snd wsIghU ^ sptisiisl 



UNCLASSIFIED 



1473 (BACK) 



39 



Security Classificgtioa 



DISTRIBUTION LIST 



Documents Department 
General Library 
University of California 
Berkeley, California 9^720 

Lockheed-Callfomia Company 
Central Library 

Dept. 77-14, Bldg. 170, Pit. B-1 
Burbank, California 91503 

Naval Ordnance Test Station 
China Lake, California 
Attn: Technical Library 

Serials Dept., Library 
University of California, San Diego 
La Jolla, California 92038 

Aircraft Division 

Douglas Aircraft Ccmpany, Inc. 

3855 Lakewood Boulevard 
Long Beach, California 9O8OI 
Attn: Technical Library 

Librarian 

Goveriment Publications Rocm 
University of California 
Los Angeles 2 k, California 

Librarian 

Numerical Analysis Research 
University of California 
405 Bilgard Avenue 
Los Angeles 24, California 

Chief Scientist 
Office of Naval Research 
Branch Office 
1030 East Green Street 
Pasadena, California 9IIOI 

Commanding Officer and Director 
U. S. Navy Electronics Lab. (library) 
San Diego 52, California 



General Dynamlcs/Convair 
P. 0. Box 1950 
San Diego, California 92112 
Attn: Engineering Ilbrsury 

Mail Zone 6-I57 

Ryan Aeronautical Company 
Attn: Technical Information 

Services 
Ilndber^ Field 
San Diego, California 92112 

General Electric Company 
Technical Information Center 
P. 0. Drawer QQ 

Santa Barbara, California 93102 
Library 

Boulder Laboratories - 
National Bureau of Standards 
Boulder, Colorado 

Government Documents Division 
Uudversity of Colorado Libraries 
Boulder, Colorado 80304 

The Library 

United Aircraft Corporation 
4 00 Main Street 

East Hartford, Connecticut 06108 

Documents Division 
Yale University Library 
New Haven, Connecticut 

Librarian 

Bureau of Naval Weapons 
Washington 25, D. C. 

George Washington University Library 
2023 G Street, N. W. 

Washington, D, C, 20006 

National Bureau of Standards Library 
Room 301, Northwest Building 
Washington, D. C. 20234 



40 



Director 

Naval Research Laboratory 
Washington, D. C. 20390 
Attn: Code 2027 

University of Chicago library 
Serial Records Department 
Chicago, lUlnols 60637 

Documents Department 
Northwestern University Library 
Evanston, Illinois 

The Technological Institute, Library 
Northwestern University 
Evanston, Illinois 

Librarian 
Purdue University 
Lafayette, Indiana 

Johns Hopkins University Library 

Baltimore 

Maryland 21218 

Martin Ccmpany 
Science-Technology Library 
Mail 398 

Baltimore, Maryland 21203 

Scientific and Technlceil Information 
Facility 

Attn: NASA Representative 

P. 0. Box 5700 
Bethesda, Maryland 20014 

Documents Office • 

University of Maryland Library 
College Phrk, Maryland 207^2 

The Johns Hopkins University 
Applied Physics Laboratory 
Silver Spring, Maryland 
Attn* Document Ubrarlan 

LLbrarian 

Technical Library, Code 24-5L 
Building 39/3 
Boston Naval Shipyard 
Boston 29) Massachusetts 



Massachusetts Institute of Technology 
Serials and Documents 
Hayden Library 
Cambridge 39) Massachusetts 

Technical Report Collection 

3O3A, Pierce Hall 

Harvard University 

Cambridge 38, Massachusetts 

Attn: Mr. John A. Harrison, Librarian 

Alumni Memorial Library 
Lowell Technological Institute 
Lowell, Massachusetts 

librarian 

University of Michigan 
Ann Arbor, Michigan 

Gifts and Exchange Division 
Walter Library 
University of Minnesota 
Minneapolis, Minnesota 55455 

Reference Department 
John M. Olln Library 
Washington University 
6600 MLllbrook Boulevard 
St. Louis, Missouri 63130 

Librarian 

Forrestal Research Center 
Princeton University 
Princeton, New Jersey 08540 

Engineering library 
Plant 25 

Grumman Aircraft Engineering Coip. 
Bethpage, L'. I. , New York Il7l4 

Librarian 

Fordham Ifalverslty 
Bronx 58, New York 

U, S. Naval Applied Science Laboratory 
Technical library 
Building 291, Code 9832 
Naval Base 

Brooklyn, New York H25I 



41 



Librarian 

Cornell Aeronautical Laboratory 
4455 Genesee Street 
Buffalo 21 , New York 

Central Serial Record Dept. 

Cornell University Library 
Ithaca, New York l4850 

Columbia University Libraries 
Documents Acquisitions 
555 W. H4 Street 
New York, New York 10027 

Ehgineering Societies Library 

545 East 47th Street 

New York, New York IOOI7 

library-Serials Department 
Rensselaer Polytechnic Institute 
Troy, New York 12l8l 

Librarian 

Documents Division 
Duke University 
Durham, North Carolina 

Ohio State University libraries 
Serial Division 
1858 Nell Avenue 
Columbus 10, Ohio 

Commander 

Rilladelphla Naval Shipyard 
IMladelphla, Pennsylvania I9II2 
Attnt librarian, Code 249c 

Steam Engineering Library 
Westinghouse Electric Corporation 
Lester Branch Post Office 
Philadelphia, Pennsylvania 19115 

Hunt Library 

Carnegie Institute of Technology 
Pittsburg^ 15, Pennsylvania 

Documents Division 
Brown University library 
Providence, Rhode Island 02912 



Central Research Library 
Oak RLdge National Laboratory 
Post Office Box X 
Oak Ridge, Tennessee 

Documents Division 
The Library 

Texas A & M University 
College Station, Texas 77845 

librarian 

LTV Vought Aeronautics Division 
P. 0. Box 5907 
Dallas, Texas 75222 

Gifts and Exchange Section 
Periodicals Department 
University of Utah libraries 
Salt lake City, Utah 84U2 

Defense Documentation Center (DDC) 
Cameron Station 
Alexandria, Virginia 22514 
Attn: IRS (lO copies) 



FOREIGN COUNTRIES 

Exchange Section 
National Lending library for 
Science and Technology 
Boston Spa 
Yorkshire, England 

The Librarian 

Patent Office library 

25 Southampton Buildings 

Chancery Lane 

London W. C. 2., England 

librarian 

National Inst, of Oceanography 
Woimley, Oodalmlng 
Surrey, England 

Dr. H. Tiger schiilld. Director 
library 

Chalmers University of Technology 
Oibraltargatan 5 
Gothenburg S, Sweden 



42 



Mr. John M. Goodman 
Search Radar Branch 
Radar Division 
Naval Research Laboratory 
Washington, D. C. 20390 

Dr. G. P. Ohman 
High Resolution Branch 
Radar Division 
Naval Research Laboratory 
Washington, D. C. 20390 

Mr. Allan Oster 

Space Research Div. (Code N223.3) 
Astronautics Dept. 

Naval Missile Center 
Point Mugu, California 

Mr. Jerald Leish 

Avionics Division (Code 3222) 

Systems Dept. 

Naval Missile Center 
Pt. Mugu, California 

Mr. Edward Keller 
Bureau of Naval Weapons 
Code (RTOS 311) 

Washington, D. C. 

Dr. A. B. Dember 
Astronautics Dept. (Code N22) 
Naval Missile Center 
Point Mugu, California 

Mr. W. A. Eberspacher 

Space Research Div. (Code N223) 

Astronautics Dept. (18 copies) 

Naval Misfeile Center 
Point Mugu, California 

Mr. D. Sullivan 

Technical Director (Code NOl) 

Naval Missile Center 
Point Mugu, California 



Capt. C. O. Holmquist 
Commander, Naval Missile Center 
(Code NOO) 

Point Mugu, California 

Mr. A. Solferino 
Space Systems Dev. Div. 

(Code N222) Astronautics Dept. 
Naval Missile Center 
Point Mugu, California 

Dr. Leonard Porcello 
Radar Laboratory 
Institute of Science & Techn. 
University of Michigan 
Ipsilanti, Michigan 

Dr. William Brown 
Radar Laboratory 
Institute of Science & Techn. 
University of Michigan 
Ipsilanti, Michigan 

Mr. S. Weisbrod 

Smyth Research Association 

3555 Aero Court 

San Diego, California 92123 

Mr. J. Smyth 

Smyth Research Association 

3555 Aero Court 

San Diego, California 92123 

Mr. Robert Anderson 
Space Research Div. (Code N223.4) 
Naval Missile Center 
Point Mugu, California 

Mr. William Titus 
Avionics Div. (Code 3222) 

Systems Dept. 

Naval Missile Center 
Point Mugu, California 



43 



Li 



Mr. Edward Orn stein 
High Resolution Branch 
Radar Division 
Naval Research Laboratory 
Washington, D. C. 20390 

Dr. Irving Page 
Radar Division 
Naval Research Laboratory 
Washington, D. C. 20390 

Mr. A. Conley 

Bureau of Naval Weapons 

Code (RTOS 41) 

Washington, D. C. 

CDR J. R. High 

Director, Astronautics (Code N2) 
Naval Misfeile Center 
Point Mugu, California 

Chief of Naval Operations 
Codes OP-76, OP-321 (2 copies) 
Washington, D. C. 

Chief, Bureau of Naval Weapons 
Code RT 

Washington, D. C. 



16 t Lti 66 



0 I SPLAY 



TA7 
.U62 
no. 56 



TA7 
.U62 
.no. 56 



Chaney 



8341 



Electron ballistics 
and electromagnetic 
waves in the ionosphere. 



r LJ 



DISPLAY 



83412 

Chaney 

Electron ballistics 
and electromagnetic 
waves in the ionosphere. 



7.U62no.56 
Electron ballistics 







